Geometric Phase in Entangled Bipartite Systems 
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The geometric phase (GP) for bipartite systems in transverse external magnetic fields is investi- 
gated in this paper. Two different situations have been studied. We first consider two non-interacting 
particles. The results show that because of entanglement, the geometric phase is very different from 
that of the non-entangled case. When the initial state is a Werner state, the geometric phase is, 
in general, zero and moreover the singularity of the geometric phase may appear with a proper 
evolution time. We next study the geometric phase when intra-couplings appear and choose Werner 
states as the initial states to entail this discussion. The results show that unlike our first case, the 
absolute value of the GP is not zero, and attains its maximum when the rescaled coupling constant 
J is less than 1. The effect of inhomogeneity of the magnetic field is also discussed. 

PACS numbers: 03.65. Vf, 03.65.Ud 



I. INTRODUCTION 



A system can retain the information of its motion when 
it undergoes a cyclic evolution, in the form of a geometric 
phase (GP), which was first put forward by Pancharat- 
nam in optics [1| and later studied explicitly by Berry 
in a general quantal system @|. Great progress has been 
made in this novel region The original adiabatic con- 
dition in Berry's work has been removed by Aharonov 
and Anandan [4( , and Samuel and Bhandari have general- 
ized the geometric phase by extending to noncyclic evolu- 
tion and sequential measurements |5| • At the same time a 
kinematic approach to the theory of the geometric phase 
has also been developed by Mukunda and Simon [fj . Re- 
cently the generalization to mixed states was conducted 
first by Uhlmann in the mathematical context of purifi- 
cation |jj , and then by Sioqvist et al. based on the Mach- 
Zender interferometer |8| . Consequently the mixed-state 
geometric phase has been experimentally verified using 
both NMR interferometry [j| and single photon inter- 
ferometry [lOj. Recently the geometric phase for mixed 
state was put forward by Singh et al. for the case of de- 
generate density operators [ll| and a general formula for 
the parallel transporting was also provided. Despite the 
great progress in this field, Bhandari recently raised the 
criticism that the definition for geometric phase in mixed 
state fails when the interference fringes disappear |12| . 
This can be explained as the disappearance (appearance) 
of the geometric phase (off-diagonal geometric phase). 
The definition of off-diagonal geometric phase(OP) was 
first given by Manini et al. for pure states in adiabatic 
evolutio ns 1131 , and then was generalized to non-adiabatic 
situation[lJ| and in mixed-states|15| . Moreover the off- 
diagonal geometric phase was studied in degenerate case 
^(|and in bipartite systems 01 • Recently the effect of en- 
tanglement on the off-diagonal geometric phase has also 
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discussed |l8| . 

The quantum computation scheme for the geometric 
phase has been proposed based on the Abelian 0] or 
non-Abelian geometric phase |2dj |. in which geometric 
phase has been shown to be intrinsic against faults in the 
presence of some kind of external noise due to the geo- 
metric nature of Berry phase. Consequently the quantum 
gates based on geometric phase have also been proposed 
in different systems [2l|. where the interactions play an 
important role for the realization of some specific opera- 
tions. 

Bipartite system is of great importance in quantum 
computation, such as the transfer of quantum informa- 
tion, the construction of entanglement as well as the re- 
alizations of logic operations. Furthermore, it was found 
that GPs may be used to design quantum logic gates. 
These facts together give rise to a question what are the 
geometric phase and its motion in bipartite systems. Re- 
cently some papers have addressed this issue l23j . 
where the discussions respectively focus on the entan- 
glement dependence of the geometric phase for subsys- 
tem and the coupling effect on the GP for subsystem 
under adiabatical evolution. However, another question 
remains open; how the entanglement or interaction af- 
fects the geometric phase for the whole system. This 
consideration is not trivial since any quantum informa- 
tion procession cannot implemented by only one qubit. 
Moreover some interesting results may be found with the 
increment of the size of the system. In this paper, focus- 
ing on the entanglement and inter-subsystem couplings, 
we present explicit discussion on the geometric phases in 
bipartite systems. 

For this purpose, our discussion is divided into several 
sections. In Sec. II, we describe the model to be dis- 
cussed in this paper and some formulas are also present 
for the calculation of the geometric phase. Then in Sec. 
Ill we study the geometric phase for two noninteracting 
particles and some interesting results can be found in this 
section. The intra-subsystem coupling effect on the geo- 
metric phase is studied in Sec. IV. Finally, conclusions 
are given at the end of this paper. 
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II. MODEL 

In order to highlight the effect of the entanglement, we 
choose the initial state as 

p(0) = i^/ + r|$)($|, (1) 

where r £ (0, 1] determines the mixing of this state and / 
is the unit matrix in the 2x2 Hilbert space. Notice that 
for r = Eq. JIJ is the unit matrix and its geometric 
phase is trivial (always be zero), our discussion excludes 
this case. The state |4>) may be either of the following 
states, 

\if) =sm0|ll)i,2 + cos0|OO)i, 2 

\j>) =sm0|lO)i,2 + cos0|Ol)i, 2 (2) 

where 9 determines the degree of entanglement and 
|l(0))i(i = 1,2) is the eigenstate of Pauli operator a z . 
One should note that when 9 = (3)7r/4, the equations 
above are Bell states and Eq.Q) are the Werner states 
|24j , which plays an important role in quantum informa- 
tion processings, especially in quantum communication 
via noisy channels ,[25j and quantum distillation [2q . In 
general the mixed state Eq. which first was intro- 
duced by Wootters [2j} , may be entangled: that can be 
judged by the Peres-Horodecki condition |2g. Eq. (JIJ in- 
cludes all possible cases, such as pure or mixed state and 
maximal or non-maximal entangled states. One should 
note that Eq. is triplet-degenerate for r ^ 1. 

We should point out that the initial state Eq. (|TJ is not 
a trivial generalizing from the pure case. The first term in 
Eq. (JIJ can be regarded as the noise and the mixing coef- 
ficient r properly describe the intensity of noise. Recently 
the one-to-one correspondence between r of Werner state 
and the temperature T of the one-dimensional Heisen- 
berg two-spin chain with a mag netic field B along the 
z-axis, has been established 29]. This connection gives 
us the strong physical support for the initial state Eq. 
©• 

We choose the system composed of two spin-1/2 par- 
ticles undergoing spin procession in an external time in- 
dependent magnetic field in the z direction. Then the 
Hamiltonian is 

H = H + Hj, (3) 

in which the free Hamiltonian Hq and the XX interaction 
Hj are respectively 

Hq = LU1SI+OJ2S2 

Hi = giStSz+S^Sl), (4) 

where g > is the antiferromagnetism coupling constant 
and Sf(i = 1,2) is the z component of spin operator re- 
spectively. Sf — Sf + iSf{i = 1,2) are the raising and 
decreasing operators of the z component of spin-1/2. Ac- 
tually the Hamiltonian Eq. © is two-qubit XX model, 



which is of fundamental importance to understand the 
relation between the entanglement and quantum corre- 
lation in interacting many-body systems. In general we 
suppose that u)\ may not be equal to L02 because of the 
inhomogeneity of the external magnetic field. One will 
find in the following discussion that the inhomogeneity of 
the external magnetic field has nontrivial effect on GP. 

Some formulas should be addressed for the calculation 
of GP in this model. Recently GP for the mixed state 
has been discussed by Sjoqvist et al.in based on the 
Mach-Zender interferometer and a formula was provided 
for calculation of GP for a mixed state, 

7 9 = ArgTr[C/ll(i) p (0)], (5) 

111 which [/H(t) was defined as the parallel transportation 
operator. Based on this work, Singh, ef al., |ll| studied 
GP for non-degenerate and degenerate mixed states and 
provided a general method for the construction of [/"(f) 
by imposing 

U\t) = U{t)V{t) (6) 

in which [/(f) is the unitary evolution operator and is 
equal to e~ lHt in our model and V(t) is a blocked matrix, 
of which the elements are determined by 

\fj), \v), |//), \v') S {degenerate subspace} 
Vkk — e*(*'' ff ' fe ) t , \k) £ {the remaining space}, (7) 

where |/z), \v), |//), \k) are the eigenstates of p(0), 
and the interference terms between the degenerate space 
and the other space are set to be zero in order to keep the 
parallel transporting in the degenerate space. One notes 
that for different mixed states one has different C/H(f). 
Based on these formulas, one can calculate GP. In the 
following calculations, we label the initial states as 

1 - r 

Pi = -^-h + r\(p)(ip\ 

P2 - + »#)(# (8) 

and the geometric phases for these states are respectively 
calculated in the following parts. 

III. g = CASE 

We first study the geometric phase without interaction. 
It should point out that this situation still has interest in 
quantum information, such as the quantum teleportation 
and communication in which the nonlocal correlation, i. 
e. the entanglement between two space-liked particles 
plays a crucial role. The geometric phases for pi(i=l,2)in 
this case can be easily obtained from the formulas in the 
former section, 
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-r(sin[cos 29^^-t] cos m&Ht - cos 29 cosfcos sin ^4^t) 
,y — arctan L 1 ± _ ^ ' ± ' 

' 9 ±== + I^( CO s[cos26>^f^]cos^iii^^ 



-rfsin cos 20^=^ cos ^-^i - cos 2(9 cos cos 29^^-t] sin 
7q2 = arctan^ \— 2 i i i 2 J 2 > / g v 

+ ^ cos cos26»^=^t cos^ a t + cos2^sincos26lii^^2t sini i;i f^t) 




FIG. 1: The geometric phases 7 9 j(j = 1,2) [Arc] versus r, 9. 
For y g i (a), we have chosen u>i = ui2 and io\t = n/2; whereas 
for 7 9 2 (b), loi — 2ui2 and uit = 7r. 
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FIG. 2: The geometric phase [Arc] vs the inhomogeneity of 
magnetic fiedl n = with r = 1. Since the figures is sym- 

metrical with the 9 — ir/2, we only draw for 9 € [0, 7r/2]. (a) 
corresponds to pi and the other parameters are the same as 
that of Fig. Q[a). The dotted, dashed, longer-dashed and solid 
lines correspond respectively to n — 0.99, 0.9, 0.5, 0. When 
n — 1, the geometric phase is not discontinued at the point 
9 — 7r/4. (b) is for p2 and the other parameters are the same 
as that of Fig. H^b). The dotted, dashed and solid lines corre- 
spond to n — 0.5, 0.1, 0.01. When n = 0, the geometric phase 
is not discontinued in the point 9 = tt/4. (c) and (d) demon- 
strate the case that there is no the singularity ( We have 
chosen ui\t = 7r/4 for pi and ui\t = it/2 for p2 respectively). 



It is interesting to note that when the initial states are 
Werner states, the geometric phase is zero when the de- 
nominators in Eqs.© are not vanishing simultaneously. 
Furthermore when uj\ = u>2, the geometric phase for pi 
is zero since pi is commutative with the Hamiltonian 
and it cannot pick up any geometric phase. A detailed 
demonstration for jgjij = 1,2) with the parameters r, 9 



is shown in Fig. ^ From the figures it is obvious that 
the absolute values of 7 s i(2) attain the maximum when 
6 i= 7r/4, 3tt/4 and because of the mixing, scaled by r, 
the absolute value of geometric phase is compressed and 
tends to be zero with r — * 0. We also note that a singu- 
larity about "f g i appears when 9 = (3)vt/4 and r = 1, as 
displayed in Fig. [fl[a). At this point, the numerator and 
the dominator in the expression of j g i in Eq. Q are si- 
multaneously zero and the geometric phase is undefined. 
One has to calculate the so-called off-diagonal geomet- 
ric phase to retain the information of the evolution |13| . 
Furthermore, our calculation shows that the singularity 
depends not only on the entanglement, but on how the 
system evolves. 



With the consideration of the two noninteracting par- 
ticles, it is of interest to discuss the effect of the inho- 
mogeneity of the external fields. The results have been 
illustrated in Fig. [21 It clearly shows that when there 
is a singular point, for p%, the absolute value of the ge- 
ometric phase with n — > 1 attains the maximum closed 
to this point. Whereas, for p2, this happens for n —> 0. 
We should emphasize that this phenomenon appears only 
when there exists the singularity, which is induced by 
the entanglement of the initial state. If there is no pos- 
sibility of the appearance of the singularity, our calcu- 
lation shows that the points where the absolute value of 
the geometric phase attains the maximum is independent 
of the inhomogeneity of the external fieldjust shown as 
Fig|2Ic)and (d) and always is zero for the Werner state. 
Physically it is thus possible for this phenomenon to act 
as the signature of the singularity of the geometric phase, 
which is usually originated from the degenerace |13|. 



Another interesting consideration is that one supposes 
U)2 = 0, which corresponds to case that particle A is pro- 
cessing with frequency u>i while particle B keep dynami- 
cally fixed. In this case, the affect from the dynamics of 
the second particle is excluded and one can more clearly 
check the effect of the entanglement on the geometric 
phase. The geometric phase for this special situation can 
be obtained by setting w 2 = in Eq. ©. One can find 
that because of the entanglement, the geometric phase 
for particle 1 is zero too for Werner state, and it is also 
possible for the singularity to appear, for example when 
r = 1 and uj\t = tt. 
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FIG. 3: GP for 7 9 2[Arc] versus r and rescaled coupling con- 
stant J. We have chosen 9 = 7r/4 (a) and 9 — 3n/4 (b) for 
the highlight of the couplings effect on the geometric phase, 
which corresponds to the Werner state as the initial state, 
and set u)\t = 2u 2 t = ir. 



IV. g CASE 

Because of the intra-subsystem coupling, the evolu- 
tion of the system can be very different from the free 
case. So in this section we focus on the effect of coupling 
on the geometric phase. Obviously one note that since 
[Hi, pi] — 0, Hi has trivial effect on 7 9 i and thus this 
state have been excluded in this section. Based on the 
formulas in Sec. II, GP are for p 2 , 
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arctan 



"(sin X\t cos X 2 t — cos X\t sin X\t) 



l-r 
2 



l+r 

2 



(cos Xit cos X 2 t + sin X\t sin Xit 



in which 



Ai = [(wi - w 2 )cos26» - 25sin26»]/2 
A 2 = VK-^2) 2 +V/2. 



The results have been illustrated 
rescaled coupling constant J = — 



in Fig. with 
Compared with 
the g = case, we have chosen the Werner states 
(6 = (3)-7r/4) as the initial states and edit — 2uj 2 t = tt 
with consideration of the inhomogeneity of external mag- 
netic field. Because of the coupling, the geometric phase 
for the Werner state is obviously very different from the 
free case. From these figures we also note that GP has 
a maximal or mininmal values when the rescaled cou- 
pling constant is less than 1. Furthermore, when J > 1 
the absolute value of GP decrease with the increment of 
J. One also notes that because of the mixing, which is 
scaled by the parameter r, the absolute values of GP is 
compressed, which is same as the free case. The effect of 
the inhomogeneity of magnetic field also discussed. Our 
calculation shows that with lo 2 /lui — * 1, the point that 
GP has maximal or minimal value is infinitely closed to 
J = 0.5 (Fig. EJ. 

We also investigate the GP when one particle in zero 
field. This case is nontrivial since any system cannot 
avoid the affect from the other party. Also it can be used 
to manipulate the behavior of one particle by changing 
of the coupling strength. To this end, we set u 2 = 
and the geometric phase have been illustrated in Fig. [5] 
Different from Fig. [3| the absolute value of GP attains 
the maximal value when J — > 1. 
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FIG. 4: GP[Arc] vs. the inhomogeneity of magnetic field 
n — W2/W1 . We have choose the Werner state (9 = 7r/4)as 
the initial state with r = l,uit = 7r.The inset is for n = 0.8. 
Similar behavior can be found when 9 — 37r/4. 
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FIG. 5: The geometric phase [Arcjfor p 2 when u>2 = versus 
r, J. Except of u>2 = 0, all other settings are the same as in 
Fig. |3 



The interesting extension to this discussion above is in 
the limit of large J = — . In this case GP has a novel 
character, which have been displayed in Fig. H3 We note 
that GP tends to be zero in the limit of large J. This 
can be understood as that in the limit of large J the 
interaction Hi is dominant in Eq. 10 and the Hamil- 
tonian is inclined to be commutative with the Werner 
state. Then the entanglement of the initial state tends 
to be unchanged. 

Based on the analysis above, we can conclude that 
because of the inter-subsystem couplings, the geometric 
phase shows two different characters; in the weak cou- 
pling limit, the interaction obviously benefits the geo- 
metric phase, displayed in Fig. and Fig. However, 
with further increment of the couplings, the GP decreases 
and tends to be zero infinitely. This phenomena can be 
explained easily only if one note that in the infinite cou- 
pling limit, Hi is dominant and its eigenstates are just 
two of the Bell states. The entanglement is destroyed in 
the weak couplings limit and then revived by the interac- 
tion in the infinite limit. The interesting relation between 
entanglement and interaction between two parties is an 
important aspect of quantum information. 
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J J 

FIG. 6: 7 S 2[Arc] in the limit of large J. The parameters have 
same values to that of Fig. |3]and r = 1 is chosen. 



V. CONCLUSIONS 

In conclusion, we have discussed the geometric phase 
for entangled mixed state Eq. in an external mag- 
netic field. For the free case (g = 0), our studies show 
that because of the entanglement, the geometric phase 
for the system displays two different types of character. 
The first is that if there was no singularity, the absolute 
value of the geometric phase attains the maximum when 
the initial state is not the Werner state, independent of 
the inhomogeneity of the external field (see Fig. J2Jc) 
and (d)). Because of the mixing, scaled by r, it is com- 
pressed. Furthermore the geometric phase is always zero 
for the Werner state. The second behavior occurs when 



the singularity appears, induced by the entanglement of 
the initial states under the proper evolution. In this case 
the geometric phase reaches the maximum with 9 tend- 
ing to the singular point (see Fig. da) and Fig. [2Ja) 
and (b)). We also discuss the situation that one particle 
is processing and another keeps dynamically fixed, and 
similar results can be found. 

For the case g ^ 0, we choose the Werner states as 
the initial condition for facilitating our discussion. The 
results show that GP is completely determined by the 
rescaled coupling constant J and attain the maximal or 
minimal value when J < 1. Similar to the free case, we 
also discussed the effect of the inhomogeneity of the ex- 
ternal magnetic field. The results show that for Werner 
state GP is always zero when system is in a homogeneous 
external magnetic field , independent of the interaction. 
Furthermore with u^/wi — > 1, the absolute value of GP 
attains a maximum at the point J — > 0.5. The geomet- 
ric phase for LU2 — is also studied. Further study in 
the limit of large J displays a novel phenomena that GP 
tends to be zero and this result have been explained as 
the revival of entanglement in this limit. 

This work was supported by NSF of China under grant 
10305002 and 60578014. 
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